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Abstract 



We compute the non-Gaussianity of the curvature perturbation generated by ekpy- 
rotic collapse with multiple fields. The transition from the multi- field scaling solution 
^ . to a single-field dominated regime converts initial isocurvature field perturbations to an 

almost scale-invariant comoving curvature perturbation. In the specific model of two 
fields, (j>i and 4>2, with exponential potentials, — V{ exp(— Cj0j), we calculate the bispec- 
trum of the resulting curvature perturbation. We find that the non-Gaussianity is dom- 



£S) \ inated by non-linear evolution on super-Hubble scales and hence is of the local form. 

The non-linear parameter of the curvature perturbation is given by f^L = — 5c|/12, 
where c,- is the exponent of the potential for the field which becomes sub-dominant 
at late times. Since cj must be large, in order to generate an almost scale invariant 
£ — , | spectrum, the non-Gaussianity is inevitably large. By combining the present observa- 

tional constraints on /nl and the scalar spectral index, the specific model studied in 
this paper is thus ruled out by current observational data. 

X 

1 Introduction 

The existence of an almost scale-invariant spectrum of primordial curvature perturbations 
on large scales, with an approximately Gaussian distribution, is one of the most important 
observations that any model of the early universe should explain. The inflationary scenario 
offers a possible explanation in terms of vacuum fluctuations of scalar fields during an ac- 
celerated expansion preceding the standard hot big bang, though attempts to embed such a 
scenario within a fundamental theory such as superstring/M-theory may require some degree 
of fine-tuning (see PQ for a recent review). 
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The ekpyrotic scenario [2] (see also [31 H] ) is one alternative approach where the large-scale 
perturbations are generated from vacuum fluctuations during a collapse phase driven by a 
scalar field with a steep, negative exponential potential. It was shown, however, that even 
though the Bardeen potential acquires a scale-invariant spectrum during ekpyrotic collapse, 
the comoving curvature perturbation has a steep blue spectrum in the original ekpyrotic 
scenario driven by a single scalar field [5] . If the contracting pre-big bang phase is connected 
to the expanding hot big bang phase through a regular four- dimensional bounce, then we 
expect the comoving curvature perturbation to remain constant for adiabatic perturbations 
on large scales [6], [8] , which means that the growing mode of curvature perturbations in 
the expanding phase also acquires a steep blue spectrum. 

One way to avoid this is to consider non-adiabatic perturbations, which require two 
or more fields 0, [TO]. Recently, there has been progress in generating a scale-invariant 
spectrum for curvature perturbations in the ekpyrotic scenario with more than one field, 
which we will refer to as the new ekpyrotic scenario [TTJ [T21 [13] . If these fields have steep 
negative exponential potentials, there exists a scaling solution where the energy densities 
of the fields grow at the same rate during the collapse [HI [15]. In this multi-field scaling 
solution background, the isocurvature field perturbations have an almost scale-invariant 
spectrum [13], owing to a tachyonic instability in the isocurvature field. The multi-field 
scaling solution in the new ekpyrotic scenario can be shown to be an unstable saddle point 
in the phase space and the stable late-time attractor is the old ekpyrotic collapse dominated 
by a single field [T6] . 

The existence of a tachyonic instability raises questions about initial conditions in the 
new ekpyrotic scenario [17], but the transition from the multi-field scaling solution to the 
single-field-dominated solution also provides a mechanism to automatically convert the initial 
isocurvature field perturbations about the multi-field scaling solution into comoving curva- 
ture perturbations about the late-time attractor [18]. In this case, the final amplitude of the 
comoving curvature perturbation is determined by the Hubble scale at the transition and if 
this parameter and the initial conditions are set appropriately, the prediction for the primor- 
dial curvature perturbations from this scenario is compatible with an almost scale-invariant 
spectrum (see [TTl [T2"| [T3"l [19] for other mechanisms to convert the initial scale-invariant 
isocurvature perturbations into curvature perturbations). 

Recently, the non-Gaussianity of the distribution of primordial curvature perturbations 
in the inflationary scenario has been extensively studied by many authors (see e.g. [2U] for 
a review). Measurements of the non-Gaussianity already provide important constraints on 
specific models of the early universe and such measurements will continue to improve in the 
near future, for instance with the Planck satellite [21]. Since the non-Gaussian signal from 
single-field, slow-roll inflation is suppressed by slow-roll parameters to an undetectable level 
|22j . if non-Gaussianities are detected then this would rule out many models of inflation. 

Thus, as a natural extension of the study performed in [TBI HE] , in this paper we compute 
the non-Gaussianity of the primordial curvature perturbations generated from the contract- 
ing phase of the multi-field new ekpyrotic cosmology. We adopt the same specific model 
as the previous analysis of Refs. [161 EH] where the scale-invariant curvature perturbation is 
generated by the transition from a scaling solution, with two fields driving the collapse, to 
a single-field dominated regime. We assume that the ekpyrotic collapse is subsequently con- 
verted to expansion by a regular bounce, during which the comoving curvature perturbation 
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is conserved on large scales and the curvature perturbation generated during the collapse is 
thus directly related to the amplitude of the observed primordial density perturbation. 

This paper is organized as follows. In Sec. [2] we briefly review the background dynamics. 
In Sec. [3] we summarize the statistical properties of linear and non-linear perturbations. The 
<5iV-formalism is introduced which is used to compute the primordial curvature perturbation. 
In Sec. H] we study linear perturbations while in Sec. [5] we generalize this study to non- 
linear perturbations and quantify the expected non-Gaussianity. In Sec. E] we draw our 
conclusions. In three appendices we review previous results for the linear fluctuations of 
the iso curvature field during the ekpyrotic phase and calculate its intrinsic non-Gaussianity 
using the interaction Hamiltonian, as well as presenting a numerical check of our analytical 
results. 



2 Homogeneous dynamics 

We first review the background dynamics of the fields in the new ekpyrotic cosmology with 
multiple scalar fields. During the ekpyrotic collapse the contraction of the universe is assumed 
to be described by a 4D Friedmann equation in the Einstein frame with n scalar fields with 
negative exponential potentials 

3H 2 = V + jr , (1) 

3 

where 

n 

V = -J2 Vje-^i , (2) 

3 

and we take Vi > and set SnG equal to unity. 

From now on, for simplicity, we concentrate our attention on the case of two fields. In 
this case, it will be easier to work in terms of new variables [16J, 



C 2 4>1 + Ci^2 Ci0i - C 2 02 fo . 

<P= — / = X= — / = ( 3 ) 



cj + cl \Jcl + cl 

corresponding to a fixed rotation in field space. The potential given by Eq. (j2J) can then be 
simply re-written as jTH [161 123] 

V = -U( X )e- c r (4) 

where 

C 3 C i 

and the potential for the orthogonal field is given by 

U( x ) = Vi e~ (ci/c2)c * + V 2 e (c2/ci)cx , (6) 

which has a minimum at 



c\ + cl \<%V 2 , 



X = Xo = ^^ln(i-^) . (7) 
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If we expand U(x) in Eq. ([6]) about its minimum we obtain 



u( x ) = u 



(? CC^ 



where 



ci + 4 



(9) 



The multi-field scaling solution corresponds to the classical solution along this minimum 
X = Xo> while ip is rolling down the exponential potential. The explicit form of the multi-field 
scaling solution is given as 



(-ty, 

^lnH)-ilnf^M' 
c c \ U 



(10) 
(11) 



where p = Y,j 2/ c ? = 2/c 2 . The potential for x has a negative mass-squared around x — Xoj 



d 2 V 
dx 2 



c 2 V < 



(12) 



and thus x represents the instability direction. Furthermore, the x field evolution is nonlin- 
ear, with the cubic interaction being given by 



d 3 V 
dx 3 



cm 2 

A 



(13) 



which becomes important when we consider the non-Gaussianity later in this paper. 

If the initial condition for x is slightly different from xo or x is not zero, then x starts 
rolling down the potential and the solution approaches a single-field-dominated scaling so- 
lution. The explicit form of the single-field-dominated scaling solution is given as 



-t) 

— ln(-t) ln' Fn ■ 



(14) 
(15) 



where pj = "2/c 2 . In this paper, we consider the case in which the background evolves from 
the multi-field scaling solution to the 02-dominated scaling solution without loss of generality. 



3 Statistical correlators 

Here we briefly summarize the statistical properties of the scalar field fluctuations during 
ekpyrotic collapse. Then, in order to link these to the observable primordial curvature 
perturbation we discuss the 5iV-formalism. 

In the two-field new ekpyrotic cosmology, the isocurvature fluctuations acquired by the 
field x during the multi-field scaling regime, play a crucial role to generate a scale-invariant 
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spectrum of perturbations. On the other hand, the fluctuations of the field <p are negligible 
on large scales, because of its very blue spectral tilt [HI [TJl [T3] . Thus, in the following we 
neglect Sip fluctuations. 

In linear perturbation theory, when interactions are neglected, the free-field fluctuations 
Sx are Gaussian. Their power spectrum V x can be defined by 

2vr 2 

(5 Xkl 5 Xk2 ) = (2 7 r) 3 5 3 (k 1 + k 2 )^-V x (h) , (16) 

where the angle brackets denote an ensemble average. 

If the cubic self-interaction in Eq. (8) is taken into account Sx is no longer Gaussian and 
the first signal of non-Gaussianity comes from the three-point correlation function. Similarly 
to Eq. (fTBT) . the bispectrum of Sx, B x , is defined by 

(<WXk 2 <W = {2n) 3 6 3 Q2%) B x( k uh,h). (17) 

j 

To characterize the bispectrum one can also define the nonlinear parameter of the field 
fluctuation, f^ L , as 

When the intrinsic non-Gaussianity is local in real space this parameter is k independent 
and Sx can be written as 

8x = 8xi. + lftM> (19) 

where 5xl is the linear and Gaussian part of the field fluctuations. 

To relate the non-Gaussianity of the scalar field fluctuations to observations, we need 
to calculate the three-point functions of the comoving curvature perturbation (. In order 
to do that, we can use the 5iV-formalism [231 129 123 123 I2HI I22j- In the (^-formalism, 
the comoving curvature perturbation ( evaluated at some time t = tf coincides with the 
perturbed expansion integrated from an initial flat hypersurface at t = £j, to a final uniform 
density hypersurface at t — tf, with respect to the background expansion, i.e., 

C(t f , x) ~ 6N(t f , U, x) = N(t f , t h x) - N(t f , U) , (20) 

with ^ ^ 

Af(t f ,U,x) = 1 1 H(x,t)dt, N(t f ,U) = f f H(t)dt, (21) 

where TC(x., t) is the inhomogeneous Hubble expansion. We can calculate SN on large scales 
using the homogeneous equations of motion, in the assumption that the local expansion on 
sufficiently large scales behaves like a locally homogeneous and isotropic universe, according 
to the so-called "separate universe" approach [l3D3 ED E] • This allows us to compute the full 
nonlinear curvature perturbation in the large-scale limit. Note that we leave the initial time 
ti unspecified and we are free to identifying it with the time of Hubble crossing t — i*, i.e. 
the time when a mode k exits the Hubble radius during inflation, k = aH, or with a later 
time. 
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We will chose the initial time t{ to be during the multi-field scaling regime. Furthermore, 
since cp is unperturbed, 5N can be expanded in series of the initial field fluctuations 5\i- 
Retaining only terms up to second order, we obtain 

5N = N^5 Xl + l -N^ x% {5 Xl )\ (22) 

where N jX denotes the derivative of N with respect to x- 

Now we can convert the higher-order information about the initial field fluctuations into 
the statistical properties of the observed primordial curvature perturbations. The power 
spectrum of the comoving curvature perturbation (, V^, is defined as 



(C kl Ck 2 ) = (27r) 3 5( 3 )(k 1 + k 2 )^-V c (h) . (23) 



At lowest order, from Eqs. (|T6l) . ([22"]) and (1231 . is expressed as 

n = Nf Xi V Xi . (24) 

Note that N yXi is independent of wavenumber k and hence the scale dependence of the 
primordial spectrum, An = dhiV^/dliak, is given by the spectral tilt of the field fluctuations, 
An x , on the initial hypersurface [given in Eq. (169]) in Appendix IA"] . 

The bispectrum of the curvature perturbation (, which includes the first signal of non- 
Gaussianity, is defined as 

(C kl Ck 2 Ck 3 > = (2n) 3 6^C£^)Bdh,k 2 ,k 3 ), (25) 

3 

where the left hand side of Eq. ( |25l) can be evaluated by the 5iV-formalism using Wick's 
theorem, 

(CkiCkaCks) = N a i ( S Xik 1 SXik2 5 Xik 3 ) + l N f Xi N ,xiXi( S Xik^Xik 2 (SXi * ^Xi)k 3 ) + perms . (26) 

In the above equation, a star * denotes the convolution and we have neglected correlators 
higher than the four-point. 

Observational limits on the non-Gaussianity of the primordial curvature perturbations 
are usually given on the nonlinear parameter f^L defined by [22] 

6 f B < (27) 

If the non-Gaussianity is local, one can write ( as 

8N = CL + lfNL&, (28) 
5 

where (l is a Gaussian variable. 
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To compute the bispectrum of the curvature perturbation one can use the <5iV-formalism 
and after some manipulations, from Eqs. (1251) and (1261) one finds 



k N 

B c (k x , k 2 , k 3 ) = N^B^h, k 2 , k 3 ) + 4vr 4 P c 2 ^^ , (29) 

where B Xi (ki, k 2 , k 3 ) is the bispectrum of the scalar field evaluated at t = tj. Together with 
Eq. fT27"|) the nonlinear parameter f NL becomes 

f _ /ivx i 5 N jXiXi 

'As ,Xi 

where we have used Eqs. ffl9|) and fT27"j) . The first term on the right hand side of Eq. (1301) . 

fXi 

A3) _ JjVL /oi\ 
J NL — ) \ OL ) 

comes from the intrinsic three-point correlation functions of the field fluctuations and con- 
tains also the non-Gaussianity of quantum origin generated on small-scales, i.e., inside the 
Hubble radius during the ekpyrotic collapse described by the multi-field scaling solutions. 
The second term on the right hand side of Eq. (1301) , 



f (4) _ 5 N ^ /0O\ 

fNL ~ Q N 2 . W 

:Xi 



is completely momentum independent and local in real space, because it is due to the evo- 
lution of nonlinearities outside the Hubble radius during the multi-field ekpyrotic collapse. 
Note that the splitting in fffi L and fffi L depends on the time tj. 



4 Linear curvature perturbation 

The power spectrum of x is given by (see [IB] and Eq. (170)) in Appendix A) 



(33) 



in the limit of large e, where the fast-roll parameter e is defined as e = —H/H 2 = c 2 /2 = 1/p, 
and it is constant for the multi- field scaling solution. Thus, in this limit, V x has a scale 
invariant spectrum. 

On the other hand, as shown in [TB] in the fast-roll limit and assuming an instantaneous 
transition from the multi-field scaling solution to the single-field 02- dominated scaling solu- 
tion, the power spectrum of the final (after the transition) curvature perturbations becomes 



4 + 4 



Hi 



2n 



(34) 
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where Hx is the Hubble parameter evaluated at the transition time t = tx- The scalar 
spectral index of ( is 

An = A[\ + \) . (35) 



Note that it is always positive, i.e., the power spectrum is always blue 

Let us now interpret the result above in the light of the <5iV-formalism. We consider 
the situation in which x% is perturbed on the t = ti hypersurface, while Hi assumes on this 
hypersurface a constant value. This is justified by the fact that the t = ti hypersurface is flat 
and since x is an isocurvature field its fluctuations do not affect the local Hubble expansion. 
Furthermore, we assume that the transition into the single-field-dominated scaling solution 
at the time t = tx, happens instantaneously on the hypersurface X = Xt = const., where Ht 
is perturbed. 

Under these assumptions, the expansion N defined by Eq. ( T2~Ti) can be split into 

N = f T Hdt+ f 1 Hdt, (36) 

where tf is set sufficiently later than the transition time tx- In Eq. fl36l) . the first integral is 
over the multi-field scaling evolution and the last integral is over the </> 2 -dominated phase. 

Since the multi-field scaling solution is characterized by Eq. ffTUl) . the first term on the 
right hand side of Eq. (136|) can be expressed as (1/e) Yn.(Hi/ Hx)- Similarly, since the single- 
field dominated scaling solution is characterized by Eq. ffl4|) . the second term becomes 
(1 / €2) Itl(Hx / H f) , where e 2 = c|/2. Then, for a fixed ti and tf, the expansion N can be 
expressed as 

2 

iV = — 5 In lifrl + const., (37) 

Cl 

which depends only on the parameter C\, besides the transition time tx- 

During the multi-field scaling regime, the linear evolution equation of x on large scales 
is given by 

X + SH X + m 2 x x = 0, (38) 

where the mass of x is defined in Eq. f JT2|) . During the multi-field scaling the evolution of 
X is dominated by the tachyonic mass and thus 3Hx «C m^x- With = —2/t 2 (see the 
appendix, Sec. one finds x oc 1/t oc if, and thus 



H T = H^. (39) 

Xi 



Using this relation one can derive iV with respect to Xi an d obtain 

2 2 H T 

4xi 4xt h, 



= 7277 = 3— IT- ( 4 0) 



In particular, by using this equation and comparing Eq. (|24J) with Eq. (J34j), one obtains the 
value of xt, i-e., 



2Jci + c 2 2 

XT = V 2 2 - (41) 



In appendix C, we have calculated N Xi numerically and checked Eq. (|41|) 
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5 Non-Gaussianities 



In this section we compute the non-Gaussianity of the primordial curvature perturbation 
generated by the nonlinear dynamics during multi-field ekpyrotic collapse. We will use the 
<5iV-formalism to calculate SN(x) in the instantaneous transition approximation. We will 
thus consider quadratic terms in 5\i in the <5iV-expansion ( 1221) and, in contrast with Sec. HI we 
also allow for non-linear evolution of x on large scales. For another example of the multi-field 
calculation of non-Gaussianity from fields with exponential potential see [32], even though 
in the context of assisted inflation. 

Including the cubic self-interaction given in Eq. (TlBl) . the large scale evolution equa- 
tion for x in the multi-field scaling regime becomes 

X + 3Hx + m 2 x x = -^cm 2 x x 2 - (42) 

The above evolution equation can be solved perturbatively. Given the solution to the linear 
equation fl38l) . i.e., xl oc H, the growing- mode solution for x is 

X = Xl + \cxl = aH + ^5a 2 H 2 , (43) 

where a is a constant parameter whose value distinguishes the different trajectories. Note 
that this is a perturbative result, i.e., it is valid only as long as cx <C 1. However, since x 
grows during the collapse, unless prevented by a bouncing phase, eventually this condition 
is violated. 

For perturbations in the value of x on a hypersurface of uniform H we have 

S X = (l + ^caH^j H5a + ^cH 2 (6a) 2 . (44) 

The non-linear self-interaction of x i n Eq. (41) grows in time on super-Hubble scales, and 
thus the intrinsic non-Gaussianity of x increases. Therefore, unless we take t; sufficiently 
early that we can neglect these nonlinearities, we cannot naively use Eq. (39), which is only 
valid at linear order, and its derivative with respect to Xii to estimate the non-Gaussianity 
of the curvature perturbation. We need to work with a variable that is as close as possible 
to a Gaussian. It turns out that it is convenient to choose a as such a variable. Indeed, if we 
assume that 5xl, the perturbation of the solution of the linear equation (|38|) . is Gaussian, 
then also 5a is a Gaussian random variable because 5a = 5xl/H. Comparing this equation 
with Eq. (fl9|) we find that the intrinsic non-Gaussianity of x is °f local form and time 
independent, and it is given by 

f X NL = ^c. (45) 

In Sec. [B] of the appendix, we have checked that this result agrees with the one obtained 
using the approach of Maldacena [22] with an interaction Hamiltonian containing [13J, 
and confirms that we can consistently assume 5xl, and thus 5a, to be Gaussian. The reason 
for this is that the classical nonlinear evolution on large scales completely dominates over 
the sub-Hubble nonlinear interactions of quantum nature. 
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We will now compute the non-Gaussianity of the curvature perturbation using Eqs. (13"0T 
. The final result for f NL will be manifestly independent of the time we choose to evaluate 
5xi although the splitting into f^' L and is dependent on tj. 

Since 5 a can be assumed to be Gaussian, the simplest way to compute f^L is to calculate 
the 5N corresponding to the fluctuation 5a, i.e., 

5N = N, a 5a + ^N, aa (5a) 2 . (46) 

In order to compute iV jQ , and N j0ca we want to use Eq. (1371) . and for this we need to know 
how Ht varies as a function of a at the transition from multi-field scaling to single-field 
02-dominated scaling solution. Inverting Eq. (H3|) (to leading order in c\) gives 

<. = !(!- 'fx) . (47) 

Assuming as in the linear case that the transition corresponds to a critical value of the 
tachyon field \ = xt, on the transition surface (constant \t) we have from ( 14T1) that a oc Hrp 
and hence we find 

r , r 2 fo 1 

5N = - ^ — , 48 

cf a cf \a J 

which means 

2 1 9 1 

N = — - N = -— — (49) 
cf a cf a 1 

Taking 5a to be a Gaussian random variable and comparing with Eq. (128]) with Cl = 
—25a/(c\a) we obtain the nonlinear parameter for the curvature perturbation after the 
transition: 

5 N aa 5 2 /r- n \ 

fNL= 6J^ = -12 C ^ (50) 

This is our main result. The non-Gaussianity is given in terms of ci, where — V\ exp(— ci<pi) is 
the potential of the field 0i which remains subdominant after the transition from the multi- 
field scaling to the single-field dominated regime. The non-Gaussianity is of local form. This 
is due to the fact that it is generated by the nonlinear super-Hubble evolution. 

Equation fl5T?|) includes also the non-linear growth of the tachyon field on large scales due 
to its self-interaction. In order to see this, we can compute f§' L and fff L defined in Eqs. (1311) 



and (1321) . Thus, we have to identify the linear and non-linear dependence of 5N on the field 
values 5\i on the initial hypersurface. 

Replacing Eq. (T4"5l) in Eq. fl3Tl) we obtain 

f ^ L= u^ t ' (51) 

Furthermore, from Eqs. ( 1471) and ( 14H1) we have 

dN da 2 / 1_ \ 

N » -^^i 1 - ~H • (52) 
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and hence 

/S = ^c X , (53) 

- U*-*>¥ T < (54) 

where for the last equality we have replaced c using its definition, Eq. (jUJ), expressed x% i n 
terms of xt using the linear relation Eq. fl39l) . and replace xt using Eq. (jHJ). 
Secondly we have the contribution due to 

d 2 N ( daV dNd 2 a 2 

' x *"^teJ + ^^!"~^f' (55) 

Note that this relation is valid to linear order in cXi, and thus contains also the nonlinear 
self- interaction of x generating non-Gaussianities after t{. Using Eq. ( 1321) . this gives 



5 2 ^ / 2 2\ -"i /C7\ 

= cf Cn — c; — . 57 

12 1 12 v 2 l7 F T v ; 

Both fffi L and depend upon the choice of the initial hypersurface and thus of \%- However 
the total f^L comes from the sum of the two terms, it is independent of ti, and is given by 
Eq. (J5DJ. 

In appendix C, we calculated N >a and N >aa numerically and confirm that the result 
Eq. (130]) obtained by the instantaneous transition and the fast roll approximations is satisfied 
with good accuracy. 



6 Conclusion 

In this paper we have studied the nonlinear evolution of perturbations in the multi-field new 
ekpyrotic cosmology. If one sets the model parameters and initial conditions appropriately, 
then the prediction for the power spectrum of curvature perturbations produced in multi- 
field ekpyrotic collapse can be constrained by present observations [HI Q2J [131 EE] • In order 
to distinguish this model from other early universe scenarios, such as inflation, by future 
observations, it is important to estimate the non-Gaussianity of the curvature perturbations. 

We have studied the simplest model based on two fields with exponential potentials 
and considered the specific scenario in which the nearly scale-invariant comoving curvature 
perturbation is generated by the transition from the multi-field scaling solution to the single- 
field dominated attractor solution. We have applied the <5iV-formalism, which is widely 
adopted to study the non-linearity of the primordial curvature perturbation from inflation 
[28] . to ekpyrotic cosmology. We identify the non- linear curvature perturbation on uniform- 
density hypersurfaces at late times with the perturbed local expansion, 5N, with respect 
to an initial spatially flat hypersurface. The primordial non-Gaussianity parameter f^i 
is a sum of two contributions: f^\i defined in Eq. ( 131 p . comes from the intrinsic three- 
point function of the isocurvature field perturbation 5x on an initial hypersurface, and /Jvx, 
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defined in Eq. ( 132|) . originates from the nonlinear relation between the primordial curvature 
perturbation and the isocurvature field perturbations during multi-field ekpyrotic collapse. 
It should be emphasized that although the decomposition of Jnl into f$ L and f^ L may be 
convenient, it is unphysical and depends upon the choice of the initial time tj. However, we 
show that the physical quantity, the total f^L, is independent of t{. 

Both the multi-field and single-field ekpyrotic solutions are power-law solutions. We find 
a general result that in case of a sudden transition between two power-law solutions the local 
expansion is only a function of the Hubble rate at the transition, H?- Thus the calculation 
of the primordial curvature perturbation in our model reduces to finding the perturbation 
of Ht for different trajectories in phase-space, and hence the local value of the isocurvature 
field x on the initial spatially flat hypersurface. 

We find that after the transition to the single-field attractor solution the non-Gaussian 
parameter f NL = — 5c^/12, where — V\ exp(— Ci0i) is the potential of the field cf>i which 
becomes subdominant at late time. Since the non-Gaussianity is mainly generated by the 
nonlinear super-Hubble evolution, it is of the local form, and the nonlinear parameter is k 
independent. We show in appendix [B] that the contribution of the intrinsic non-Gaussianity 
of the isocurvature field perturbations on sub-Hubble scales is subdominant. 

We have checked our analytical results by calculating the expansion, N, numerically in 
appendix We confirmed that when the fast-roll parameter satisfies 400 > e > 25 the 
analytic estimation is accurate within at a few % level. The discrepancy arises from the 
breakdown of both the fast-roll approximation and the sudden transition approximation. 

A negative value of /nl is much more tightly constrained by current observations than 
a positive value. For instance, if we choose C\ = 5 we obtain a nonlinear parameter /nl ~ 
—10 that is marginally consistent with current constraints on the non-Gaussianity of the 
primordial density perturbation [33] [M]. However such a small value of c\ leads to too large 
a value of the spectral index (35), An > +0.16, which is excluded by observations [33J. 

Corrections to the exact exponential potentials that we have studied in this paper have 
been proposed [TJJ [T2] to produce a red spectrum of perturbations. Indeed corrections are 
also required to stabilise the multi-field scaling solution at early times (before observables 
scales exit the Hubble scale) [17]. We assume that the classical background solution starts 
close to the multi-field scaling solution. However, as the multi-field scaling solution is a saddle 
point in the phase space, we need some preceding phase that initially drives the classical 
background solution to the saddle point. It is important to check whether such corrections 
will also affect the non-linear evolution of the field perturbations during the collapse phase, 
which could lead to additional sources of primordial non-Gaussianity. 

Modifications to the effective action are certainly required at high energies to turn con- 
traction to expansion before the collapse phase reaches the big crunch singularity [351 ttB 
[T2"l [T3] . In this paper, we assumed that the comoving curvature perturbation is conserved 
on super-Hubble scales through such a non-singular bounce, neglecting the effect of non- 
adiabatic perturbations which are expected to rapidly decay about the single-field attractor 
solution. In this case the power spectrum and the non-Gaussian parameter Eq. (150]) that we 
have calculated in the single-field dominated collapse are directly related to those observed 
perturbations in the expanding, hot big bang phase. It is certainly possible to convert the 
initial isocurvature field perturbations to curvature perturbations through a different mech- 
anism, e.g., via the bounce [11] [12] [13] [T7J. In this case the non-linear dependence of the 
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expansion, 8N, upon the initial field perturbations may be different, but the non-linear 
growth on super-Hubble scales of the isocurvature field perturbations about the multi-field 
scaling solution is still expected to lead to a large non-Gaussianity [13j El- Further work is 
required to quantify the non-Gaussianity predicted in these models. 

Note added: This arXiv version of our article includes corrections made after publication 
of the journal article. We give the correct expressions for several equations that have an 
incorrect sign. This incorrect sign is due to a sign error in the integrated expansion N 
in Eq. (37), which then propagates into other equations. As a consequence, the nonlinear 
parameter /nl in this new ekpyrotic scenario turns out to be negative and it is much more 
tightly constrained by current data than it would have been with a positive sign. By combining 
the present observational constraints on / NL and the scalar spectral index, the specific model 
studied in this paper is thus ruled out by current observational data. 
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A Linear field fluctuations 

Here, we briefly summarise the results obtained in the previous works pXjQjZJQ^Q^EIE] about 
the generation and the evolution of linear fluctuations of the x field about the background 
multi-field scaling solution. 

In the multi-field scaling solution 5\ coincides with entropy field perturbation, which is 
automatically gauge-invariant and does not couple with gravity at first order. Since there is 
no coupling between the adiabatic and the entropy field perturbations, neglecting the non- 
linear self-interactions of the entropy field which are negligible at sufficiently early times, 
the equation of motion for 5\ becomes simply the equation of a free massive field in an 
unperturbed FRW metric, 



where m 2 has been defined in Eq. ( 1X21) . Introducing the rescaled field v = a5x, and writing 
the wave equation in terms of conformal time r = J dt/a, we have 




(58) 



v, TT + k 2 - 



a 



i 2 2 

+ m v a v 







(59) 



a 



where (. . .) iT denotes the derivative with respect to r. 
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For the multi-field scaling solutions, we can show that the following relations hold, 

aH = — — — , (60) 
(e-l)r 

^n = _( e _2)a 2 tf 2 , (61) 
a 

m 2 = 2e(3 - e)H 2 . (62) 



Thus, Eq. ( T58|) becomes 



2 , e + 3r/ x - 2 
(e - l) 2 r 2 



A; 2 + 



u = , (63) 



where r\ x = m 2 /(3if 2 ). 

Using the usual Bunch-Davies vacuum state to normalise the amplitude of the fluctuations 
on small scales, we obtain 

- i(i/+l/2)f 

^-kT^H^i-kr), (64) 



2 A; 1 / 2 

where the order of the Hankel function is given by 

At late times, — kr — > 0, making use of the asymptotic form of the Hankel function, 

Hi>H-k T )^-m(=^y , (66) 

the power spectrum of 8\ in this limit becomes 

V x = cA-kr) 1 ^ , (67) 
where C v = 2 z/ ~ 3 / 2 r(z/)/7r 3 / 2 . Then the spectral tilt of the generated fluctuations is 

and to leading order in a fast-roll expansion, 

2 

An x ~ - . (69) 

Thus, for a steep exponential potential we obtain a slightly blue spectrum, becoming scale- 
invariant as e — > oo. 

Approximating ^ ~ 3/2, on large scales we can relate the amplitude of b% to H as 

= ^ • (70) 
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B Isocurvature field bispectrum 

In this section we calculate the intrinsic non-Gaussianity of the field Xi using the approach of 
Maldacena [22] , which also includes the contribution from sub-Hubble nonlinear interactions. 
Following [22], the three-point correlator is given by 

(5 X 3 (t)) = -i f dt([5 X 3 (t),H int (t>)]) , (71) 

J — oo 

where H int is the interaction Hamiltonian. Here we consider only the cubic interaction, which 
is the dominant one for the three-point function, 

r V® 

H int (f} = J d 3 xa 3 ^-5x 3 , (72) 

where is given in Eq. (fl3|) . 

Writing Eq. (ITT]) in Fourier space, 

3 

ft 



j 

2Re (-i<W*)<W*)<W*) f ^ dt'a 3 V® 6x1^x1^x1^')) > ( 73 ) 



6 Xk = -^=e- ikT ( 1 - ± ) (71) 



and using the normalised free field solution 

1 

a\f7k V kr 
valid in the limit p <C 1, we get 

y( 3 ) e 2 h-4 

(*x kl (*)W*)<W*)) = -^m^k,)— ^ 

3 J J 



r f ^ (1 - zAur') (1 - ik 2 r') (1 - ik 3 r') k > T '] , (75) 

J — oo T 



i (1 + ifor) (1 + ik 2 r) (1 + i/c 3 r) e 
dr' 

OO 7 



where we have taken out from the time integral using that /H 2 is constant. 

Taking into account only the leading contribution in — kjT <C 1, the last two lines of the 
above equation yield Y.jk 3 /A. Furthermore, using also Eqs. ( TT3l) and (1621) to replace V^ 3 \ 
we finally find 

<W0*XiJ*)<W*)> = (2ir) 3 6(£k j )^^H i . (76) 

3 3 J 

The intrinsic non-Gaussianity after Hubble-exit is thus of the local form, as it is dominated 
by the super-Hubble evolution. Note that using e = 1/p and H = —p/t we can rewrite the 
last term in the equation above as 

- w < 77 > 
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Thus, for c = 1/M we recover the result of [13]. With the definition of the power spectrum 
for 5x, Eq. dMD, we can rewrite Eq. (ITS"]) as 

(<W*)<W*)<W*)> = ^(E^f^vr 4 ^, (78) 



which yields, using Eq. ([I 



5x(^^ 2 ^3) = |^27r 4 ^. (79) 

This equation is equivalent to Eq. (1451) . confirming the result found from the large scale 
nonlinear evolution in Sec. Thus non-linearities on sub-Hubble scales are completely 
subdominant with respect to the nonlinear classical super-Hubble evolution. 



C Numerical results 

In this appendix, we check the validity of the analysis of Sees. @] and by numerically 
calculating the expansion 5N. 

C.l Phase space variables 

For numerical calculations, it is more convenient to adopt phase space variables. Thus, first 
we summarise the background dynamics in terms of the phase space variables. These are 
defined as [13 [361 EZ] 

*> = jb> < 80 > 

* = V^T' (81) 

and their evolution equations are given by 

% = -3*i(l-E4)-^yf/|, (82) 

|| = Vif3E^-^i) > (83) 
where N = log a and j = 1,2 for two fields case. The Friedmann equation gives a constraint, 

E-'-j-E/Zj '• (84) 
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Assuming J2j c j 2 < 1/6, there are 4 fixed points of the system where dxj/dN = dyj/dN = 0, 



£4 = 1. y, = 0. (85) 



c /c 2 

■Bj : Xj = -j=, y j = -^-±-1, Xk = y k = 0, (for j ^ k) , (86) 



where p = £)j 2/c 2 . A linearized stability analysis shows that the multi-field scaling solution, 
B, always has one unstable mode. On the other hand, the single-field-dominated scaling 
solutions, Bj, are always stable. 

In the (xi,x 2 ) plane, the three fixed points B,Bi and B2 are connected by a straight 
line, which is given by 

ClC2 fQQ\ 

c 2 xi + cix 2 = —f= ■ (88) 
V6 

Solutions which start close to the saddle point B evolve along this line to one of the single-field 
dominated solutions, B\ or B 2 . As in the main text, we concentrate on the case in which the 
background evolves from the multi-field scaling solution B to the single-field 02-dominated 
scaling solution B 2 for simplicity. 



C.2 Numerical scheme 

We first briefly summarise the numerical scheme. For the phase space variables Xj, jjj defined 
by Eqs. (I8"U|) and flSTl) . we solve the evolution equations (1821) and numerically. Since 
we concentrate on the case in which the background evolves from the multi-field scaling 
solution, initial conditions are characterised by a point close to B in the phase space. As in 
the main text, in the multi-field scaling solution background, the fluctuations of the field ip 
are negligible. 

Thus, for the phase space variables, we choose the initial conditions 

Vqh 

-x 

V6H 

where Eq. (189]) is given by the background scaling solution in Eqs. ( 1101) and (fTTj) . This 
guarantees that (^1,^2) lies on the line given by Eq. (1881) . 

Equation (|90|) defines the quantity z, which characterises the deviation along the insta- 
bility direction from B. Once we fix a model, i.e. the values of C\ and c 2 , the initial values 
of Xi and x 2 are determined from Eqs. (|89p and (|9"U|) for a given initial value of z. 

For this set of (xi,X2), yi is fixed so that (x±, X2,Vi) is in the instability direction from 
B, and then y 2 is fixed from the constraint equation (184"|) . Using this set of initial values 
(xi,x 2 ,yi,y 2 ), we can specify a background trajectory by solving the evolution equation. 



C 2 Xi + CiX 2 



4 + 4 



c 



C 2 X 2 — C\X\ 



4 + 4 



(89) 
(90) 
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The left panel of Fig. [T] shows an example of the background solution. For this background, 
we evaluate an expansion 5N from the hypersurface characterised by Hi to that by Hf. 

From Eqs. (|43|) . (jUJ) and (|90|) . we can see that a different choice of z corresponds to that 
of a and \ as 

z = ^=(l + ~caH)a, (91) 

= + (92) 

Then Z{ is related with Xi an d it determines the Hubble parameter at the transition Ht in 
this scheme. In order to calculate the derivatives of N(zi) with respect to Zi, we change Zi 
slightly and calculate the difference of the expansion N(z{). The right panel of Fig. [1] shows 
the effect of the change in Z{ on the Hubble parameter. 




0.01004 0.01006 0.01008 



Figure 1: Left: An example of background solution for log \ H\ with initial condition z — Zi. 
We also show log|i?r| which is determined from the amplitude of 8\- Right: The same 
background solution shown around N ~ 0.01005 which is much later than the transition. 
We also show the solutions with slightly different initial condition, z = + 5zi, z = Z\ — 5zi, 
with dashed lines. After the transition, this difference of initial z generates the curvature 
perturbations which can be evaluated as the difference of iV at H = Hf. 



C.3 Power spectrum and non-Gaussian parameter 



In the following we first check numerically the linear relation Eq. (j34p . Note that from 
Eqs. f l9Tj) and fl92|) . 5z oc 5a oc 5\ holds at the linearised level. From Eqs. ff24|) and fl33|) . the 
power spectrum of the curvature perturbation is given by 



Hi 



2tt 



(93) 



In our numerical scheme, from Eq. fl92|) we evaluate N )Xi as 



N 



y/6 



N.. 



e Njzi + dzij-Njzi-dzi) 
V6 I5zi 



(94) 
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which is valid for sufficiently small 5z{. From Eq. ( |34"|) . the power spectrum of the curvature 
perturbations can also be expressed in terms of the scalar field fluctuation 5x as 



2 



4 + 4 



h 



T 



2tt 



2 = A§i 
" 4 + 4 \Sx\l 



H 



2tt 



(95) 



As Eq. (I95p and Eq. (1931) should agree, we examine whether a quantity q defined by 

becomes close enough to 1. We find that g = 1 is satisfied within about 1 % accuracy in our 
numerical simulations. 

Next, we compute numerically /jvx. From Eq. (|9Tj) . we can write down iV Q and iV QCf in 
terms of z: 



N a = N z — = -=Ht 1 + — Zi N z , (97) 




N aa = N zz ( — ] + iV,— = ? ^ ( 1 + Zi ) AT,,, + -^JV, . ( 9i 
In our numerical scheme, we evaluate as 

AZI (fc,) 2 ' 

If the instantaneous transition approximation is valid, the non-Gaussianity parameter 
f NL is given by 

'» ^ i^f- (100) 

Thus, in order to verify the accuracy of the analytical result /tvl = — 5cf/12, we examine 
whether a quantity r defined by 

27V 

becomes close to 1 or not. 

Since we checked that the choice of 5zi does not affect the results for sufficiently small 
values, we use Szi = 10 -5 . It was also verified that the time evolution of 5\ can be described 
by that of the multi-field scaling solution up to Zi = 0.01; thus, we take Zj = 0.01. Using 
these initial conditions, we calculate f^i for various combinations of c\ and c 2 . It is worth 
noting that in terms of c\ and c 2 the fast-roll parameter is expressed as e = 44/^(4 + 4)- 
For the examples we show in Table 1, the range of the value of e varies from 25 (ci = c 2 = 10) 
to 400 (ci = c 2 = 40). The results for f^L and r are summarised in Tabled! 

From Table [1] we can see that for these parameters, the results for J'nl are accurately 
described by f^L = — 5cf/12 which agrees with the result obtained analytically based on the 
instantaneous transition approximation and the fast roll approximation in the main text. 
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Table 1: c\ and c 2 dependence of the nonlinear parameter. We also compare with the results 
obtained by the instantaneous transition and the fast-roll approximations. The deviation 
from r = 1 denotes the error of these approximations. We adopt = 0.01 and 8zi = 10~ 5 . 



Cl 


c 2 


Jnl 


r 


10 


10 


-40.6165 


0.974797 


10 


20 


-40.9285 


0.982283 


10 


40 


-41.5317 


0.996761 


15 


10 


-90.4979 


0.965311 


15 


20 


-92.6425 


0.988187 


15 


40 


-93.5294 


0.997647 


20 


10 


-161.014 


0.966082 


20 


20 


-164.976 


0.989855 


20 


40 


-169.365 


1.01619 


40 


10 


-639.909 


0.959865 


40 


20 


-659.744 


0.989617 


40 


40 


-664.996 


0.997494 



We can see that the analytical results deviate from the numerical results for smaller e 
(compare the cases with c\ = c 2 ). This is because the fast-roll approximation becomes worse 
for smaller e. We can also see that, for a fixed Cj, the analytical results again deviate from 
the numerical results for smaller c 2 . This is because the potential around B becomes flatter 
along the direction to _B 2 for smaller c 2 , and the instantaneous transition approximation 
becomes worse. 
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